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Recently, a new probabilistic "data fusion" framework based on Bayesian principles has been developed on 
JET and W7-AS. The Bayesian analysis framework folds in uncertainties and inter-dependencies in the diagnostic 
data and signal forward-models, together with prior knowledge of the state of the plasma, to yield predictions of 
internal magnetic structure. A feature of the framework, known as MINERVA (J. Svensson, A. Werner, Plasma 
Physics and Controlled Fusion 50, 085022, 2008), is the inference of magnetic flux surfaces without the use of a 
force balance model. We discuss results from a new project to develop Bayesian inversion tools that aim to (1) 
distinguish between competing equilibrium theories, which capture different physics, using the MAST spherical 
tokamak; and (2) test the predictions of MHD theory, particularly mode structure, using the H-l Heliac. 
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1. Introduction 

With the advent of large scale neutral beam heat- 
ing, magnetically confined fusion plasmas have drifted 
from the simple picture of ideal magnetohydrodynamics 
(MHD), which describes the plasma as a single, station- 
ary, isotropic Maxwellian fluid. Due principally to neu- 
tral beam heating, several tokamak experiments now boast 
plasma toroidal rotation speeds that approach the thermal 
Mach speed and have significant stored energy residing in 
the energetic particle population produced by charge ex- 
change with thermals. Motivated by these developments, a 
range of new descriptions have emerged that include ther- 
mal rotation JTJ as well as energetic particles [2j. De- 
spite this, ideal MHD is still the foundation of nearly all 
analysis. Detailed magnetic reconstruction based on this 
treatment ignores the energetic complexity of the plasma, 
and can result in model-data inconsistencies, such as ther- 
mal pressure profiles which are inconsistent with the total 
stored kinetic energy of the plasma. 

A parallel development has been the improvement in 
the diversity, accuracy and resolution of plasma diagnos- 
tics. Interpretation, however, often requires a detailed 
knowledge of the plasma equilibrium. For example, infer- 
ence of the toroidal current profile j^{t/f) from line of sight 
measurements of the polarization angle requires a knowl- 
edge of the poloidal flux if/ across the plasma. Formally, 
diagnostic forward functions relate the vector of plasma 
parameters I to the measurement vector T). For a linear 
system, such as toroidal current inference in a double null 
configuration, I and T> are normally related through a re- 
sponse matrix M with additional contributions C, such that 
T) — MI+C. Inference involves inverting this relationship 
to give plasma parameters I that are consistent with the 
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data D. A widespread technique used is least-square fit- 
ting, in which prior assumptions are included via a penalty 
term in the fit. 

The confluence of higher performance plasmas with 
diagnostic development has however led to a dichotomy: 
data is often inconsistent with ideal MHD equilibria, some- 
times misleading scientists to propose new phenomena to 
describe data "artifacts". An example is filamentation of 
flux surfaces in the Rijhnuizen Tokamak, which was in- 
ferred from fluctuations in measurements of electron tem- 
perature Q. Subsequent analysis of the Thomson scatter- 
ing detection chain revealed that with the correct photo- 
electron statistics and 2D instrument profile, similar struc- 
tures could arise from noise. |@1[5]. 

Recently, a new integrated data-modelling approach 
for inference of fusion plasma parameters has emerged 
which offers a natural framework with which to resolve ar- 
tifacts from model behaviour. In contrast to least square fit- 
ting, the Bayesian approach to inference involves the spec- 
ification of an initial prior probability distribution func- 
tion (pdf), P{I), which is then updated by taking into ac- 
count information that the measurements provide through 
the likelihood pdf P{D\I). The result is the posterior dis- 
tribution P(I\D) given by Bayes' formula 

P(I\D) = P(D\I)P(I)/P(D). (1) 

The advantage of the Bayesian approach over traditional 
inversion techniques is two-fold: (i) prior knowledge, in- 
cluding known parameter inter-dependencies is made ex- 
plicit, and (ii) as the formulation is probabilistic, random 
errors, systematic uncertainties and instrumental bias are 
integral part of the analysis rather than an afterthought. 

In this work we present initial results of forward mod- 
els of current tomography in the Mega Ampere Spheri- 



cal Tokamak (MAST), we formulate Bayesian inference 
of force balance, and we identify plans for inference of 
mode structure in the HI heliac. The paper is structured 
as follows: Sec. [2]briefiy outlines the MAST experiment 
and present first results for Bayesian inference of poloidal 
flux and current profiles from MAST motional Stark effect 
(MSE) measurements. Section [3] outlines a Bayesian in- 
ference framework for force balance, synthesizes a pres- 
sure profile posterior given the flux surfaces from MSE 
and Thomson Scattering data, and computes an estimate 
of f(ft)f'(i/f), with f(ifi) the toroidal flux function. Using 
this, we are able to quantify the impact of poloidal cur- 
rents on the equilibrium configuration. Section [4] outlines 
a Bayesian inference model for mode structure and devel- 
ops a plasma model to compute candidate modes in the H- 1 
heliac. Finally, Sec. [5]contains concluding remarks. 



2. MAST and current tomography 

The Mega Ampere Spherical Tokamak, which first 
published physics results in 2001 [6|, is one of the world's 
largest spherical tokamaks. While MAST physics and 
technology development have contributed across a broad 
range of fusion science 121 E), two properties of MAST 
help motivate this research: high performance, and pre- 
cision diagnostics. It is the combination of these proper- 
ties in a spherical tokamak, which has a relatively weak 
toroidal field and therefore large sensitivity to poloidal 
currents compared to conventional tokamaks, which helps 
motivate selection of MAST to develop Bayesian inference 
for force balance. 

Recently, both MAST neutral beam injectors have 
been upgraded to 3.8 MW. This has enabled plasma per- 
formance to be routinely lifted above f3„ « 5, which was 
reported in 2005 [9|. MAST is equipped with an array of 
precision diagnostics [8|, including a high spatial resolu- 
tion, single-time point, ruby Thomson Scattering system, a 
multi-time point Nd: YAG Thomson Scattering system with 
a coarser spatial resolution, Motional Stark effect, Charge 
Exchange Recombination, and fast magnetics IfTOl . 

Our development of Bayesian inference of the cur- 
rent profile on MAST closely follows the seminal work 
of Svensson and Werner ifTTI . In that work, the plasma 
was represented as a grid of toroidal axisymmetric current 
beams, each with rectangular cross-section and each beam 
carrying a uniform current density. In MAST, we have 
placed these beams so as to fill-out the entire limiter re- 
gion as depicted in Figure [T] The magnetic field generated 
is a summation of Biot-Savart's law over current beams. 

A key advantage to using a series of current beams 
with finite cross-section to model the plasma current (as 
opposed to a filamentary model) is that the semi-analytic 
expressions for the corresponding magnetic field and vec- 
tor potential have no singularities, even at points within the 
current beam itself. Indeed, if one were to use filaments to 
model the plasma current, there would be many singular 
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Figure 1 : Plasma beam cross-sections fill out the MAST limiter 
region in the plasma beam current model. Pick-up coils 
are shown as solid bars near the poloidal flux coils, flux 
loops as a "*" on both the centre column and poloidal 
field coils, and MSE as "*" points on a chord through the 
midplane. 



points in the calculated magnetic field within the plasma 
that would make subsequent flux-surface calculations dif- 
ficult and somewhat questionable IfTTI . 



2.1 Diagnostic Signals 

In axisymmetric devices, there exists a direct relation 
between the poloidal flux, if/, and the toroidal component 
of the current vector. Specifically, defining the cylindrical 
coordinate system (R, Z, <p) with R the major axis, Z the ver- 
tical axis and (f> the geometric toroidal angle, the following 
relation holds: 



<A(fl,Z) 



A dl = IhRAa 



where is related to the toroidal current 
Savart's law: 



(2) 



via Biot- 
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and where //o is the permeability of free-space. In Eq. (|3]l r 
is the position vector r = ReR+Ze^ipe^, which e«, e z , and 
unit vectors in the R, Z and <p directions, respectively. By 
inspection, the Ar and A z components due to j^, are zero. 
Taking the curl of Eq. (|2]i and utilizing B = V x A it can 
be shown 



B R = - 



B 7 



1 diff 
RdZ' 

RdR' 



(4) 
(5) 



where B R , B z and 7?^ are the components of the magnetic 
field. From Eqs. (|2ji-(j5]), it can be seen that if is repre- 
sented by an axisymmetric beam of rectangular cross sec- 
tion and uniform current distribution, then ifr, Br and Bz 
will vary linearly with respect to the total current going 
through that beam. The importance of this will become 
clear in the next two paragraphs. 

The main diagnostics used to infer the current distri- 
bution are pickup coils P,-, full flux loops F; and the polar- 
isation angle y, of the emitted light from neutrally excited 
species during neutral beam injection due to the motional 
Stark effect (MSE). The responses of these diagnostics to 
a current running through a beam, 7, is given by 

Pi(R, Z;I) = B R (R, Z; 7) cos(0,-) + B Z (R, Z; 7) sin(0;), (6) 

F i (R,Z;I) = y,(R,Z;l), (7) 

tan y,(R, Z; 7) = 

A B Z (R, Z;I)+AiB R (R, Z; 7) + A 2 B 4 ,(R, Z; 7) 



A 3 B Z (R, Z; 7) + A 4 B R (R, Z; 7) + A 5 BJR, Z; 7) ' 



(8) 



where 6/ is the angle between a pickup coil's normal and 
the outboard midplane, and Aq, A\, Az,A^,Ai and A5 are 
constants for the particular MSE viewing geometry. The 
notation ";" denotes that the subsequent argument, in this 
case current 7, is held constant. (By convention, the nota- 
tion "|" is reserved for probability theory.) Both P, and T 7 , 
have a linear dependence on ift, B R and B z and hence they 
are linearly dependent on the total current going through 
each beam of the current plasma beam model. The func- 
tion tan y, has a non-linear dependence on Br and B z ; how- 
ever given that both neutral beams and MSE viewing optics 
all are on the mid-plane for MAST, Ai, A2, A3, A4 are all 
approximately zero. The following approximation is hence 
used for tan 7, : 

A B z (R,Z;I)+A 1 Br(R,Z;I) 



tany(7?,Z;7) = 



A 5 73 (7?,Z;7) 



(9) 



where we have retained the term A\. Since the current 
beam plasma model does not take into account poloidal 
currents, vacuum field values for B^ are used in Eq. (|9jl to 
perform the current tomography calculations. The correc- 
tion to 7?0 due to poloidal currents is investigated in Sec. 
3D 

Thus, given that P t (R, Z; 7), F,(7?,Z;7) and 
tan ji{R, Z; 7) all have a linear dependence on the 
current flowing through 7, one may write out a generalized 



prediction vector P for all the pickup coils, flux loops and 
polarization angles, as 



P = MI + C. 



(10) 



Here, M is the response matrix of the current vector I cor- 
responding to all the plasma beams modelling the plasma, 
and C represents various other contributions to the predic- 
tion vector, which are constant relative to I. 

2.2 CAR Prior and Bayesian Inference 

In IfTTI . Svennson and Werner used a Conditional 
Auto-regressive (CAR) prior to perform Bayesian infer- 
ence for the current beam model of the plasma, and this 
is the model adopted in this paper. The advantage of this 
choice in prior is that plasma beam currents are correlated 
to each other but in a spatially localised way. Thus, this 
prior has the effect of enforcing some smoothness between 
adjacent current beams, while minimising spatially long 
ranged effects due to manipulating the current in a par- 
ticular beam. This is desirable in that such behaviour al- 
lows one to get a clear relation between individual cur- 
rent beams and diagnostic measurements. To construct this 
prior distribution, the following distribution over all cur- 
rent beams is considered 



(11) 



where the superscript T denotes the transpose, and Q is 
the precision matrix. Equation (jTTJ is proportional to a 
zero-mean normal distribution, and so the conditional dis- 
tribution of one current 7, given all others 7_,- satisfies 



pdi\I-d oc exp 



{li-YjPiihfl&n) 



(12) 



for some /Jy. Matching terms between Eq. (Hi and Eq 



12 1 shows that/3, 7 = 0, B u = -Q u /Qn and n = l/Q u . 



Moreover, it turns out the symmetry of Q implies that 
BijTj = BjiTi lfl2l . We have set all variances equal. By 
setting Bjj to correspond to 7,- having a mean that is simply 
the average of the currents of the beams horizontally and 
vertically adjacent to it, the precision matrix takes the form 



(13) 



where 1 is the identity matrix and W the adjacency ma- 
trix with Wij - 1 if i and j are horizontally or vertically 
adjacent to current element 7,, and Wy = otherwise. 

From Eq. ( 12 1, it can be seen that the CAR prior is 



constructed by exactly specifying the conditional distribu- 
tion p(7j|7_,) and subsequently transforming this expression 
into a zero-mean normal distribution for 7. CAR distribu- 
tions are discussed in detail in [ 12] and have the advantage 
that hidden parameter inter-dependencies associated with 
direct manipulation of the covariance matrix in a normal 
distribution are totally avoided in the CAR construction 
(see IfTTI for a discussion on this point). 



The likelihood distribution function P(D\T), assuming 
normally distributed errors on all measurements, is a multi- 
variate normal distribution of the diagnostic predictions P, 
given an array of diagnostic measurements D correspond- 
ing to the predictions in P. This distribution is explicitly 
represented as 



P(D|I) 



1 



(2^/ 2 |E D | 1 / 2 



exp 



1 



r(P- 



D) r £ _1 (P-D) 



(14) 



where No is the number of measurements, and £ is the co- 
variance matrix of the measurements, which is determined 
experimentally. Finally, substituting Eq. (Hi and Eq. (JT4j 
into Bayes' formula shows that the posterior distribution 
p(I|D) satisfies the following relation: 

p(I|D) oc exp [-(P - D) r £ _1 (P - D) - I r Ql] . (15) 



Recalling Eq. (IOi, it is known that P is linearly de- 



pendent on I; and thus, it can be shown that Eq. (15i 
is proportional to a multivariate normal distribution in I 
(see 1111 for details). Hence, p(I|D) has a simple analytic 
representation that be directly analysed without resorting 
to using Markov-chain Monte-Carlo (MCMC) algorithms. 
The speed afforded by the analytic nature of the posterior 
makes this analysis amenable to real-time plasma control 
applications. 

With the maximum and variance of /?(I|D) deter- 
mined, Eqs. Q and Q can be used to construct the 
poloidal flux surfaces corresponding to the contours of 
t[r{R, Z). Figure[2]shows poloidal flux surfaces from MAST 
discharge # 22254 at 320 ms using pickup coils, flux loops 
and MSE. Discharge #22254 is a deuterium plasma in a 
double-null configuration, which was heated with 3.1 MW 
of neutral beam heating and a plasma current of I p = 
800 kA. The time of 320ms is analysed here as it corre- 
sponds to a firing of the high-resolution TS system closest 
to the peak /3 for this shot. The figure shows a contour 
plot of ifr(R, Z) which is calculated from I corresponding 
to the maximum of the p(I|D) distribution. Overlaid on 
the contours are traces of the poloidal field coil cross sec- 
tions and conducting surface cross sections for the MAST 
experiment. The last closed flux surface calculated from 
the plasma beam model is outlined in bold with the corre- 
sponding EFIT last closed flux surface overlaid in purple 
for comparison. One outcome of the Bayesian approach is 
generation of pdfs from which the error to the fit can be 
inferred. Figure [3] shows the inferred safety factor, q, pro- 
file and distribution obtained by sampling the posterior 200 
times. The width of the fits is a characteristic measure of 
the error in q. 

3. Inference of force balance 

The physics goal of our work is to exploit the im- 
proved resolution of diagnostics to infer the validity of dif- 
ferent force balance descriptions. Specifically, we have in 
mind development of a framework that will validate en- 
ergetic plasma resolved force balance models. To build 




Figure 2: Poloidal flux surfaces inferred for MAST shot # 22254 
at 320 ms using pickup coils, flux loops and MSE. 



towards this goal, we commence by applying Bayesian in- 
ference to the axisymmetric Grad-Shafranov equilibrium 
description of ideal MHD force balance, J x B = VP. 

The well known Grad-Shafranov equation lfT3l can be 
written 



F(R,Z) = 



(16) 



where 



F(R,Z) = -fioRUW + hoR 2 p'M + i&MfM, (17) 



d 1 dtf/ d 2 i// 
- fM>Rm) =R--- + —. 



(18) 



For a real plasma, the presence of nonideal effects will 
mean F(R, Z) is nonzero. Our long term aim is to com- 
pute P(T\D) using Bayes formula, such that P(D\T) = 
P(D\T)P(T)/P(D). Here,the vector D comprises all the 
magnetics, motional Stark effect, Thomson scattering and 
charge exchange recombination data, while the column 
vector T contains F{R,Z) evaluated from Eq. ( 17 1 at dif- 
ferent (R, Z) across the plasma. 

An integral representation for the prior P{T) can be 
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Figure 3: Safety factor, q profile as a function of normalised 
poloidal flux found by sampling the posterior 200 times 
for shot #22254 at 320ms. The poloidal flux is nor- 
malised such that if/„ = is the magnetic axis and if/„ = 1 
is the edge. 



obtained by using the transformation property 



PCF,//(<A),/W/'(iA)x 



diUW,p'M,M)fW) 
P(M^,p'W,fWf'W) (19) 



and integrating both sides across p and f . That is 



P(T) 



-I 



dp'dffm (20) 



where the integrand dp'dff = dp' (if/)df(if/)f (iff) 
ideal MHD plasma satisfies Eq 



An 



16 1, such that the resid- 
ual force is zero, and so P{T) = <J(!F). Unfortunately, 
the joint distribution function P(j,p(if/), p'{if/), f (>//)/' (iff)) 
is non-separable, preventing direct integration of Eq. 
(20 1. In future work we will fold interdependency of 
p' and f(ip)f'(i//) together to enable calculation of 

3.1 Inference of toroidal flux function 

Some progress can be made if we assume the plasma 
obeys ideal force balance, and p'(if/) is assumed to be in- 
dependent of f(i/f)f (iff) and jfiii/f). In this instance, and 
providing we are able to estimate p'(if/) and its distribu- 
tion from measurements, then f(if/)f (if/) across the mid- 
plane can be computed through Eq. (17i. As a first step 



to inference of force balance in a real plasma, we compute 
f (iff) f (iff), integrate to find f(if/), and substitute this back 
into MINERVA. By examining the changes in the position 
of the magnetic axis and q profile from the recomputed so- 
lution, we are able to quantify the impact of poloidal cur- 
rents on current tomography without the need for a sepa- 



rate magnetic reconstruction by EFIT. This folds p'(if/) into 
j<p(iff) through the influence of poloidal currents. 
For an ideal gas, the pressure is given by 



p = n e k B T e + n,k B Ti 



(21) 



with n e , n, the electron and ion density, and T e , Tj the elec- 
tron and ion thermal temperature. At 320 ms, high reso- 
lution Thomson scattering measurements of T e and n e are 
available. The data is distributed normally about each data 
point, and a mean and standard deviation provided. To ob- 
tain a pressure estimate, we have assumed that the mean 
values of density and temperature satisfy «,/« e = 0.8, and 
Tj/T e = 1.1, as typically extracted from a charge exchange 
recombination measurement. We have also assumed the 
same distribution of data as Thomson Scattering. 

We have next computed a fit for the pressure p(if/) us- 
ing Monte Carlo simulation [ 14]. The procedure is as fol- 
lows: at each radial grid point, the inverse transformation 
method is used to generate n e , T e , n,-, Ti samples that sat- 
isfying the prescribed pdfs P(n e ), P(T e ), P(ni), P(T/). To 
compute p(if/) and its distribution we have mapped each 
sample p(r) to p(iff) using if/(r) across the inboard chord 
determined from MSE. To find a smoothed p(if/) we have 
then fitted a fourth order polynomial in normalized flux. 
The sampling process is repeated until the pdf for p(if/) no 
longer changes. 

Figure [4] shows the sampled pressure and fit for the 
pressure using N s = 2000 samples. Overlaid is the poly- 
nomial fit for each sample across the midplane. The width 
of the polynomial fit (shown in black) yields a lower case 
estimate to the standard deviation in p(if/). 




Figure 4: Inferred pressure profile p(r) across inboard midplane 
radial chord at 320 ms in # 22254. The sampled data are 
shown as points, and the black lines are fits to p(i]/) for 
each sample set using a fourth order polynomial in if/. 

Using the pressure profile fit of Figure [4] we have 
sampled p'(if/) and j$ across the midplane and computed 
f(if/)f'(if/). Figure |4] shows the mean and standard devia- 
tion for p'(tA), J0 and f(if/)f'(if/). For reference, the sched- 
uler EFIT solution for this time slice is over-plotted in bold. 
While the MSE current profile is qualitatively similar to 



EFIT, the pressure profile is very different, particularly in 
the edge and core regions This difference is principally re- 
sponsible for the difference in f(if/)f'(i[r) at the core and 
edge seen in panel (c). 

Next, we have integrated f (iff)/' (if/) to obtain f(if/) 
from panel (c). This yields a toroidal flux profile that 
varies approximately linearly in poloidal flux from f(if/) = 
-0.407 at the edge to f(if/) = -0.496 at the core. That is, 
poloidal plasma currents are paramagnetic, and increase 
Z?0 by 22% at the core. Qualitatively, this poloidal current 
should lift the on-axis safety factor « 22%, increase the 
change in poloidal flux across the plasma, but not change 
the geometry of flux surfaces or the Shafranov shift. Work 
is in progress to investigate the impact of the correction to 
f(if/) in MINERVA. 
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Figure 5: Inferred parameters across the midplane inboard radial 
chord using Monte Carlo simulation. Panel (a) shows the 
pressure gradient p'(ip), panel (b) the toroidal current and 
panel (c) the toroidal flux flux function f(tfr)f'(tfr). 



4. Inference of mode structure 

A second objective of the project is to test predic- 
tions of MHD mode theory, particularly mode structure. 
To achieve this goal we intend to develop a Bayesian in- 
ference tool to identify oscillations using the mode struc- 
ture of a collection of candidate modes. The low beta he- 
liac plasma configuration of H-1NF is ideally suited to this 
task, as there is little uncertainty in the equilibrium, and 
precision control of the helical coil currents enables access 
to a wide range of magnetic rotational transform (*) pro- 
files, as shown in Fig. [6] The combination thus offers a 
rich testing ground for Bayesian mode analysis. 

For inference of mode structure, the relevant plasma 
parameters I are the set of eigenmode-specifying parame- 
ters (frequency, mode numbers, radial structure moments) 
and D is a set of Mirnov array measurements made out- 
side the plasma. By solving the eigenmode problem in the 
plasma, and computing the vacuum-region magnetic field 
corresponding to a given plasma mode we obtain the for- 
ward model for D given I. Given some prior P(T), the 
likelihood P(D\I) can then be constructed. 

In the remainder of this section we derive the for- 



ward model for Global Alfven Eigenmodes (GAEs) in H-l 
plasmas. As a starting point to this sub-task, we compute 
GAEs in a cylindrical plasma model (coordinates (r, 6, z)) 
incorporating a vacuum region in the region r p < r < r w 
which is encased by a perfectly conducting wall at r — r w . 




Figure 6: Rotational transform as a function of average minor 
radius for the typical range of K h in H-1NF, a three field- 
period heliac with average minor radius r p = 0.2m and 
major radius R = lm. 



GAEs are discrete modes that accumulate at the mini- 
mum <d 2 . of the Alfven continuum or. = k 2 v 2 . 1151 where 

mm A || A 1 ■ 

k\\ denotes the component of the wave vector parallel to 
the equilibrium magnetic field and Va denotes the Alfven 
speed. Expanding k\\ — k ■ B and substituting into a> A = 
k 2 v 2 A yields a> A cc (n - tm) 2 /p where m,n are the poloidal 
and toroidal mode numbers respectively, and p is the mass 
density profile, as sumed to take the form p = po ( 1 - r 2 / r~\ 
with po the density on-axis. We define a> 2 ijn to be the global 
minimum of to 2 A over the interval (0, r p ), if such a mini- 
mum exists. That is, a> 2 . is the smallest 0)\(rn) with both 

' mm A v uy 

a> A (ro) = and lo 2 a (r ) > for some radius r$ such 
that < ro < r p . Figure [7] shows the dependency of 
oj m i n on a-/, for different m, n modes. We have restricted 
our attention to the frequency interval < / < 40 kHz, 
where / = a>/(2n), since we are only interested in sta- 
ble modes, and current H-l diagnostics cannot resolve fre- 
quencies above 40 kHz. Theoretically, for a given Kh, 
GAEs lie at frequencies just below the lines in Figure [7] 
accumulating at <J m ;„, with an increasing number of radial 
nodes the closer the eigenfrequency lies to u> m i„. 

We have modelled the plasma region by a stellarator 
normal-mode formulation in a cylindrical plasma [ 16l fT71 . 
corrected for non constant density profile. The eigenvalue 
equation for this problem can be written 



hp 



L0 2 M(fi, 



(22) 



where a> is given in units of inverse Alfven time r A ' = 
Bq/R ypoPoi with #o the field strength at the magnetic axis, 
and where ip is defined in terms of the radial element of the 
fluid displacement 



r^ r = exp 



i{m6 - -z) 



(23) 
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Figure 7: A plot of the minimum in the continuum a> mi „/(2n) vs 
Kjj for a selection of mode numbers (m, ri). The cross-hairs 
mark the location of the eigenmode shown in Fig. [8] 



with i = V— 1. The operators L and M are give by: 



L = - 



Id 2 d 

(n — mi) r — 

r dr dr 



(n — m t) 



1 d p(r) d p{r) m 

M = r 1 » 

r dr po dr po r l 



Ds H (n — mi) 

m 

2 



(24) 
(25) 



Equation ( 22 1 is derived by averaging over helical ripple 
and assuming large aspect ratio. The Suydam stability pa- 
rameter is given by 



D s = -^'(r)Q'(r) 
where the average field line curvature is 

Q(r) = e 2 N\ 



'{r 2 i + 2 J rtdrj 



(26) 



(27) 



Here j = y—l, 'prime' denotes radial derivatives, i = ri', 
p(r) is the normalized pressure, e is the inverse aspect ratio, 
N is the number of turns made by the helical windings and 
/?o = 2fiopo/Bg is the ratio of plasma pressure to magnetic 
pressure at magnetic axis. 

Turning now to the vacuum region, the magnetic field 
B must satisfy VxZ? = and Vfi = 0. In cylindrical geom- 
etry we can write B in terms of modified Bessel functions, 
up to a scale factor.] 18] The boundary conditions 

(28) 

(29) 
(30) 
(31) 



[[B r ]] = = 
where [[B]] = B vacuum - B plasma , together with 

6(0) = o,m*i, 
£(P) = 0,M = 1, 

B r (r h) = 0, 



relate the vacuum solution to that in the plasma. 

If a? ■ appears in the plasma and is positive, a shoot- 



ing method is used to find solutions to Eqs. (22i, (28i, 



( 3 1 1, and Eq. ( 29 1 or Eq. ( 30 1 with frequency < oj < 



a> m i„. More explicitly, we select the appropriate bound- 
ary condition for £ r at r = 0, and then 'shoot' out to the 



plasma/vacuum interface, searching frequency space in the 
region < a> < o» mI „ for eigenfrequencies that are con- 
sistent with Eq. (28 i. Figure [8] shows an example of an 
(m, ri) = (3,4) eigenmode found in this manner, located 
at the point (/<•/,, a>) indicated by the cross-hairs in Figure 
[7] The global structure seen has promising features when 
compared to recent spectral measurements of the mode: 
the mode structure is global in radial extent, and the radial 
position of the peaks nodes broadly matches observations. 





f=2.8khz 







0.10 
r[m] 



Figure 8: The radial fluid element displacement as a func- 
tion of minor radius for (m,n) = (3,4) and K h = 0.73 
at the eigenfrequency indicated This eigenmode has an- 
gular frequency to = 12.1 k rad s~' 



5. Conclusions 

In this work we have described a Bayesian inversion 
framework for inference of force-balance and mode struc- 
ture. Progress has been reported in three areas: develop- 
ment of Bayesian inversion for current profiles; inference 
of the toroidal flux profile assuming the plasma satisfied 
ideal MHD; and calculation of a Global Alfven eigen- 
mnodes in a helical-ripple averaged cylindrical plasma 
model. 

Based on working by Svensson et al. we have devel- 
oped Bayesian inversion for current profiles and poloidal 
flux contours in MAST. The model comprises a grid of 
rectangular toroidally extended current beams, and is con- 
strained by pickup coils, flux looops and MSE data. Owing 
to the placement of the MSE viewing optics the tangent 
of the polarisation angle of emitted light is linear in Br 
and Bz- Signals from the pickup coils and flux loops are 
also linear in Br and Bz- This feature means the prediction 
vector for the data is linear in the currents. We have used 
a conditional autoregresssive prior to describe the initial 
current distribution, which has the advantage of enforcing 
some smoothness between adjacent current beams while 
minimizing spatially long ranged effects. Posterior calcu- 
lations of the poloidal flux contours of MAST #22254 at 
320 ms illustrate good boundary agreement to EFIT 

In Sec. [3] we outlined an inference technique for 
the extraction of force balance, as applied to the Grad 
Shafranov equation. Using Thomson scattering data and a 



Monte Carlo simulation technique, we fitted pressure pro- 
files to the inboard midplane chord of MAST, and com- 
puted p'(if/). Assuming ideal MHD force balance we then 
calculated f(ip)f'(i/r). While the f(i/f)f'(iff) profile is simi- 
lar to EFIT, the change in f(tfr) from the core to the edge is 
22%, which is slightly larger than the EFIT value. Quali- 
tatively, this poloidal current should lift the on-axis safety 
factor by w 22%, and increase the change in poloidal flux 
across the plasma, but not change the geometry of flux 
surfaces or the Shafranov shift. At present, we are sub- 
stituting the corrected profile inferred from f(iff) into 
MINERVA to compute the effect of poloidal currents on 
the calculation of j$ and tfr surfaces. This constitutes a first 
step towards combining j^, p'(tft) and f(iff)f'(iff) to validate 
different equilibrium descriptions. 

Finally, we have developed an analysis approach for 
the inference of mode structures, and constructed a GAE 
model. The model, which is in cylindrical geometry, ac- 
counts for a non constant density profile, features helical 
ripple averaging to compute the Suydam criterion and gen- 
erate t, and includes the vacuum region. The model has 
been applied to the range of * profiles accessible to H-l 
plasmas, and the minimum in the continuum computed as 
a function of helical winding current ratio /<•/,. The eigen- 
function of a candidate (m, n) = (3,4) GAE mode was de- 
termined at an angular frequency of 12.1 krad s _1 . The 
eigenfunction has similar structure to recent observations 
in HI, including the same number of radial nodes which 
occur at similar radial positions. 

In ongoing work we are developing forward models 
of Thomson scattering and charge exchange recombination 
spectroscopy for inference of force balance in MAST plas- 
mas. We will be exploiting this information not only to im- 
prove equilibrium reconstruction, but to use Bayesian in- 
ference as a tool with which to resolve competing equilib- 
rium models. In particular, we intend to compare force bal- 
ance descriptions of ideal MHD, ideal MHD with flow (T], 
two fluid models |[T9l . and energetic fluid resolved equilib- 
ria Q. The Bayesian model offers a rigorous framework 
with which to quantify the fit and thereby elucidate the im- 
portant underlying physics. With regards mode structure, 
we intend to transform the cylindrical eigenfunction to a 
beam cross-section equilibria [20 1 and identify a candidate 
mode set which spans the different type of mode structures 
observed. The next stage of this work involves deployment 
of forward models for mode structures into MINERVA to 
compute the posterior for the mode class and parameters. 
On a longer time scale we envisage the cylindrical plasma 
model will be replaced by a fully 3D ideal MHD wave 
mode code, CAS3D [21 1. As with equilibrium modelling, 
the importance of mode analysis by Bayesian inversion is 
that it offers a rigorous framework in which to identify 
wave mode structures, some of which have deleterious ef- 
fects on plasma performance. 
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